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Abstract 

We investigate a scenario in a supersymmetric SO (10) Grand Unified Theory in 
which the fermion mass matrices are generated by renormalizable Yukawa couplings 
of the 10© 120© 126 representation of scalars. We reduce the number of parameters 
by assuming spontaneous CP violation and a 1^2 family symmetry, leading to nine 
real Yukawa coupling constants for three families. Since in the "minimal SUSY 
SO (10) GUT" an intermediate seesaw scale is ruled out and our scenario lives in 
the natural extension of this theory by the 120, we identify the vacuum expectation 
value (VEV) w R of (10, 1, 3) G T26 with the GUT scale of 2 x 10 16 GeV. In order 
to obtain sufficiently large neutrino masses, the coupling matrix of the scalar 126 is 
necessarily small and we neglect type II seesaw contributions to the light-neutrino 
mass matrix. We perform a numerical analysis of this 21-parameter scenario and 
find an excellent fit to experimentally known fermion masses and mixings. We 
discuss the properties of our numerical solution, including a consistency check for 
the VEVs of the Higgs-doublet components in the SO(10) scalar multiplets. 
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Introduction: The group SO (10) is a favourite candidate for constructing grand unified 
theories (GUTs) [I]. The special interest in such theories also stems from the fact that they 
allow for type I [2] and type II [3] seesaw mechanisms (see also [I]) for the light neutrino 
masses. Confining oneself to renormalizable SO (10) GUTs, the scalar representations 
coupling to the chiral fermion fields, which are all assembled for each family in the 16- 
dimensional irreducible representation (irrep), are determined by the relation [HIE] 



where the subscripts "S" and "AS" denote, respectively, the symmetric and antisymmetric 
parts of the tensor product. The so-called "minimal SUSY 50(10) GUT" (MSGUT) [7] 
makes use of one 10 and one 126 scalar irrep for the Yukawa couplings, to account 
for all fermion masses and mixings (8]. The MSGUT contains, in addition, one 210 
and one 126 scalar irrep [7j. This model has built-in the gauge-coupling unification of 
the minimal SUSY extension of the Standard Model (MSSM). Detailed studies of this 
minimal theory have been performed HDl HH HH H31 HH US US]; in QI1 Q3] small 
effects of the 120-plet were considered in addition. It turned out that the MSGUT works 
surprisingly well in the fermion sector, provided one neglects constraints on the overall 
scale of the light neutrino masses. This, however, proved to be crucial, since the natural 
order of the neutrino masses in GUTs is too low, namely w 2 /Mgut ~ 1-5 X 10~ 3 eV, with 
v ~ 174 GeV and the GUT scale M GUT ~2x 10 16 GeV. Thorough studies of the heavy 
scalar states [HI QH EH US [211 [22] have been used to show that this MSGUT is too 
constrained [231 121] and does not allow to enhance the neutrino mass scale to a realistic 
one [23[ I2"6] . compatible with the results of the neutrino oscillation experiments (for a 
review see, e.g., [27J). One aspect of this problem is that a seesaw scale significantly lower 
than the GUT scale spoils the gauge coupling unification of the MSSM. 

An obvious attempt to loosen the corset of the minimal theory is to add the 120-plet 
of scalars. A study in that direction has been done in [28]. Earlier works considering a 
prominent 120-plet contribution to the fermion mass matrices are found in [221 EDI ED E2] ■ 
We note that 10 © 120 alone does not give a good fit in the charged fermion sector [55] . 
Thus the 126 scalar irrep is not only needed in the neutrino sector but also for the charged 
fermion mass matrices. In that case, the mass matrices of the charged fermions and the 
neutrino Dirac-mass matrix are given, respectively, by 



The Yukawa coupling matrices H, G, F belong to the scalar irreps 10, 120, 126, respec- 
tively. The coefficients k d , K d , Kg, v d denote the vacuum expectation values (VEVs) of 
the Higgs doublet components in the respective 5O(10) scalar irreps which contribute to 
the MSSM Higgs doublet H d , the rest of the coefficients refers to H u . The light neutrino 
mass matrix is obtained as 



16 <g> 16 = (10 © 126) s © 120 AS , 



(1) 



M d 
M u 
M e 



k d H + K d G + v d F, 
k u H + k u G + v u F, 
k d H + K t G - 3v d F, 
k u H + k d G - 3v u F. 



(2) 
(3) 
(4) 
(5) 



M V = M L - M D M^M T D with M L = w L F, M R = w R F, 



(6) 
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with scalar triplet VEVs Wl and wr. The mass Lagrangian of the "light" fermions reads 



C M = -d L M d d R - u L M u u R - I L M e l R - ~v L M v [v L ) c + H.c., (7) 
with {v L ) c being the charge-conjugate of v L . 

A renormalizable 50(10) scenario: The goal of this letter is a numerical study of 
the system of 3-generation mass matrices (j2J) to ([6]), taking into account the neutrino- 
mass suppression factor f 2 /M GU T- This system does not easily lend itself to such an 
investigation because it contains many parameters, thus we use some arguments to reduce 
their number. The scenario we want to investigate is defined by the following assumptions: 

i) The Yukawa coupling matrices H, G, F are real. 

ii) We impose a 7j 2 symmetry, which sets some of the Yukawa couplings to zero and 
which is spontaneously broken by the VEVs of the 120, in particular, by k^, k u , Kg, 
Kd being non-zero. 

iii) We assume w R = M G ut, with Mgut = 2 x 10 16 GeV. 

iv) We set Wl = 0, i.e., we have pure type I seesaw mechanism. 

Let us now comment on these items. Item i) can be motivated by spontaneous CP 
violation. The Z2 of item ii) is given by 



^2 -> ~1p2 



>120 —> —^120, 



where the ipj (j = 1, 2, 3) denote the fermionic 16-plets and 0i2o is the scalar 120-plet. All 
other multiplets, not mentioned in Eq. (jHJ), transform trivially. With the Z2 symmetry of 
Eq. (jHJ), the coupling matrices have the form 



hu 






















h 3 3 



H=\ hoo , G 



( g 12 

-912 



/ fn / 13 

f 2 2 I . (9) 



V -223 J V /13 / 33 

We have used the freedom of basis choice in the 1-3 sector to set 7t 13 = 0. Of course, this 
1j2 symmetry of Eq. (jHJ) is an ad-hoc symmetry, but it enhances the importance of the 120 
because its Yukawa coupling matrix G is now responsible for mixing of the second family 
with the other twoQ Item iii) is motivated by the fact that the MSGUT does not allow 
to fix the problem of too small neutrino masses by taking wr significantly lower than the 
GUT scale [231 [2U 1251 J2H1 12H]- Thus our scenario has built in that the natural neutrino 
mass scale in the MSGUT is too low. Consequently, the neutrino mass scale has to be 



1 In Eq. ||HJ|, for the definition of the Z2 symmetry, all choices i/jj — > —ipj are equivalent. With choosing 
1P2, we anticipate the result of the fit of our scenario to the masses and mixings at the GUT scale. That fit 
gives a strong hierarchy of the elements of H, which — with Eq. ([5]) — can be formulated in the usual way 
as j /in I I/122I <S I ^33 1- Furthermore, with the convention of (jSj) it is possible to have all diagonalizing 
matrices of the charged fermion masses in the vicinity of the unit matrix. 
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enhanced by the smallness of the coupling matrix F . Item iv) is a trivial consequence 
of the previous one: for small F, type II seesaw contribution to Ai u is negligible. 

Now we tackle the problem of parameter counting. Without loss of generality, we 
assume that k<i, k u and wr are real and positive. Then we define 

H' = k d H, G'=\K d \G, F'=\v d \F. (10) 

The primed matrices have the dimension of mass. The phases of the VEVs of the 120 
and 126-plets cannot be removed. Thus we write the mass matrices as 

M d = H' + e* d G' + e iQ F', (11) 

M u = r H H' + r u e i ^G' + r F e i ^F' ) (12) 

M t = H' + r t e iil G' -3e iCd F', (13) 

M D = r H H' + r D e^ D G' -3r F e iCu F' } (14) 

M v = r R M D F'- l Ml. (15) 

The ratios r#, etc., are real by definition since we have extracted the phases from the 
VEVs. Now the counting is easily done. Since we have nine real Yukawa couplings, see 
Eq. ([9]), there are nine real parameters in H', G', F' . Furthermore, there are six phases 
and six (real) ratios of VEVs, altogether 21 real parameters. On the other hand, we have 
18 observables we want to fit: nine charged-fermion masses, three mixing angles and one 
CP phase in the CKM matrix, two neutrino mass-squared differences Am 2 tm and Arrig, 
and three lepton mixing angles. 

Suppose, we have obtained a good fit for the 18 observables. Then we still have to check 
if the fit allows for reasonable VEVs and Yukawa coupling constants. A detailed discussion 
of this issue is found in Appendix A. Here it is sufficient to note that wr = Mgut and the 
determination of r R and r F by the fit fix \v d \ and \v u \ via \v d \ — r^M GUT and \v u \ = r F \v d \. 
Therefore, first test we check 

\ Vd \ 2 + \ Vu \ 2 = \ Vd \ 2 (l + r p) < v 2 with v = 174 GeV (16) 

for every fit. Clearly, this inequality holds at the electroweak scale, and we assume that 
approximately it is valid at the GUT scale too. 

A numerical solution: To find a numerical solution, we employ the downhill simplex 
method [34J for minimizing a x 2 -function of the parameters — for an explanation of the 
method see [2T)l [55] . Actually, the x 2 -function can be minimized analytically with respect 
to the parameter r R of Eq. ffl5l) . which results in a x 2 -function depending the remaining 
20 parameters, and we apply our numerical method to that function. To build in the 
inequality ffl6|) in our search for the minimum, we add a suitable penalty function to our 
X 2 - Our scenario is fitted against the values of the 18 observables at the GUT scale; for 
an MSSM parameter tan (3 = 10, these values are displayed in Table [TJ 

Choosing the normal ordering mi < irt2 < of the neutrino masses (Am^ = m^—ml, 
Am^ tm = m| — mf), we have found a fit with af = 0.0087, which is a perfect fit for 
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Quarks 


m d 


l.DUdO _ Q 2304 


m s 


29.9454 + 4 ;300i 


m b 


1063.6 +£JJ 


m u 


0.7238 +0-1365 

— U. 14D / 


m c 


210.3273 +19-0036 


m t 


8243 3 3 +30267.6 
5Z4d^ _ 1476g 6 




0.2243 ±0.0016 


S23 


0.0351 ±0.0013 


S13 


0.0032 ±0.0005 


°CKM 


60° ± 14° 



Leptons 


m e 


n +0.0003 

U - d0S0 _o.0003 




75.6715 +0-0578 

— U.UoUl 


m T 


1292.2 + - 0013 
-0.0012 


Arn| 


(7.9 ±0.3) x 10" 5 


Am atm 


(2.2 +0-37 ) x io-3 


S 2 


0.31 ±0.025 


S 2 
*23 


0.50 ±0.065 


S 2 
tb 13 


< 0.0155 



Table 1: Input data at the GUT scale for M GU t = 2 x 10 16 GeV and tan/3 = 10. The 
charged- fermion masses are taken from [35], the remaining input from Table I in [26] . 
Charged- fermion masses are in units of MeV, neutrino mass-squared differences in eV 2 . 
We have used the abbreviations S12 = sin #12, etc. The angles in the left table refer to the 
CKM matrix, in the right table to the PMNS matrix. 

all practical purposes. This fit is so good that it does not make sense to show the pulls0 
The matrices H', G' and F' for our fit are given by 



H' 



G' 



F' 




(17) 



where all numerical values are in units of MeV; the values of the ratios of VEVs and the 
phases are shown in Table [2J The neutrino mass spectrum turns out to be hierarchical 
with m 1 = 1.57 x 10~ 3 eV < m 2 = 9.03 x 10~ 3 eV < m 3 = 46.96 x 10~ 3 eV, and the 
PMNS phasd is 12°. We want to stress, however, that our fit solution is perhaps not 



2 The largest pull is 5 x 10 2 for m s . 

3 We use the same phase convention as for the CKM matrix in 
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91.0759 


- 


- 


- 


- 




19.66974° 


r F 


297.758 


Cu 


-2.96594° 


- 


- 




189.12385° 


r u 


7.14572 


Cu 


226.65689° 


re 


1.33897 




6.24258° 


r D 


3008.88 


Co 


179.85271° 


Tr 


2.90553 x 10~ 17 







Table 2: The values of the phases and ratios appearing in the mass matrices (lllj) - (TT5]) in 
the case of our fit. Hyphens in the left two columns indicate that the ratio corresponding 
to the phase has been absorbed in one of the primed matrices, whereas hyphens in the 
right two columns signify that there is no physical phase associated with that ratio. 

unique, because with the numerical method used here we could miss other minima of % 2 . 

For our fit, it turns out that y = \va\ \Jl + r 2 F = 173.0 GeV. This looks dangerously 
close to the upper bound of Eq. ([16]) . To check if this danger is serious, we have plotted in 
Fig. [1] the minimal x 2 as a function of y. In order to pin y down to a given value y we have 
extended the x 2 function to (x 2 )y — X 2 + {(y ~~ 2/)/(0-01?/)} 2 , minimized {x 2 ) y and plotted 
X 2 at this minimum versus y — for previous uses of this method see, for instance, [26J. We 
read off from Fig. [I] that x 2 is minimal at y = 173 GeV, however, this minimum is rather 
flat; note that x 2 is plotted on a logarithmic scale. Thus we still obtain excellent fits if 
we go to lower values of y. In Appendix A, a consistency condition is worked out which 
the 5*0(10) GUT has to fulfill in order to reproduce the VEV ratios of Tabled There we 
also show that for our fit all Yukawa couplings stay in the perturbative regime. 

In order to find out if our scenario makes a prediction for the PMNS phase 5, we 
treat it in the same way as y in the previous paragraph, i.e., we consider (x 2 )<5 = X 2 + 

[(5 - 5)/ (QMS)} . Departing from 5 = 12° for our numerical solution given by Eq. (|T7j) 

and Table [21 and going stepwise down close to 5 = 0° and up to 5 = 360°, the quality of 
the fits remains excellent, with (x 2 )s always below 0.3. Thus, in our scenario all values of 
the PMNS phase are possible. 

One may ask the question how large enough neutrino masses and an atmospheric 
mixing angle which is close to maximal are accomplished with the numerical values given 
by Eq. (|T7|) and Table [21 We concentrate on achieving J Aml tm ~ 0.05 eV. With the value 
of tr we find that Tr x 10 9 MeV ~ 0.029 eV. Thus we take into account all contributions 
to M. u /tr which are of order 10 9 MeV. Rewriting Eq. (|15[) in the form 

M u = r R | (r H H' + r D e^G') F'' 1 (r H H' + r D e iio G') T - 6 r H r F e Ku H' + 9 r\ e^F'j , 

(18) 

a numerical analysis shows that all elements of the second and third term on the right- 
hand side are smaller than 1.8 x 10 s MeV and 0.5 x 10 8 MeV, respectively. Thus the 
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250 



\v d l(l+^) 1/2 [GeV] 



Figure 1: The x 2 as a function of \vj\ J 1 + t\ 



dominant matrix elements in M. v stem from the first term and are given by 



Mf°^=r R r D 



( 

r D (2g'M\, + fi 3 (g' 12 ) 
V -r H h' 33 g' 12 f[ 3 /d 



/d 



\ 

-r H h' 33 g' 12 f{ 3 /d 

r D (g 23 ) 2 h 22 



(19) 



Here, d = f'uf'^ — (/13) 2 is the determinant of the corresponding 2x2 submatrix of 
F' and we have used the approximation £d = 180°. Apart from the common factor tr, 
in this matrix there are four products of three matrix elements: one matrix element is 
always from F'~ l , the other two are either both from r^G' or one from r^G' and one from 
r H H' . Looking at Eq. (TTT1) and Tabled we find that products of the largest elements, for 
instance (g^) 2 fi,^ never occur, these would be too large. Plugging in the numerical values 
of the parameters, we find that all non-zero terms in J^i^ ora ^ are similar in magnitude: 



r R x 10 9 MeV 



-1.77 + 2.64) -0.81 
-0.81 1.00 



(20) 



Due to the minus sign in the first term, we end up with —1.77 + 2.64 = 0.87, close to 
1.00 and thus leading to nearly maximal mixing. In this crude approximation, which is 
only relevant for the largest mass and the atmospheric mixing angle, we obtain m 3 ~ 0.05 
eV and 9 23 ~ 43° Apart from the smallness of F', it is the large factor td in Md 



'Note that there is also a small contribution from Mi. 
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which gives the correct magnitude of the neutrino masses. Maximal atmospheric neutrino 
mixing rather looks like a numerical contrivance in our scenario. 

Finally, a word concerning the low-energy values of the quark masses is at order. 
Ref. [35J takes the quark mass values at mz from Ref. [37] as input for the renormalization 
group evolution up to Mgut, whereas Ref. [37] uses the input 

m u (l GeV) = 4.88 ±0.57, m d (l GeV) = 9.81 ± 0.65, m a (l GeV) = 195.4 ± 12.5 (21) 

and 

m c (m c ) = 1.302 , m b {m b ) = 4.34 +^ , m t {m t ) = 171 ± 12, (22) 

see Tables I and II in [37J. The light quark masses are given in MeV, the heavy ones 
in GeV. Comparing these values with those given in the Review of Particle Properties 
of 2006 (RPP) [36J, we find that the heavy quark masses are in reasonable agreement. 
However, in the last years the values of the light quark masses have become significantly 
lower [36] : 

m u (2 GeV) = 1.5 3.0, m d (2 GeV) = 3 7, m s (2 GeV) = 95 ± 25, (23) 

Note that one has to take into account the scaling factor mj(l GeV)/mj(2 GeV) = 1.35 
(i = u, d, s) to compare Eq. ([21]) with Eq. (123]) [36J. In order to assess the influence 
of lowering the light quark masses, we have performed a second fit, using the values of 
Eq. (1231) as input, scaled to Mgut with the factor 0.200 for m u and 0.207 for and m s 
(see [351 EZ])i but leaving the previous values for the heavy quark masses. We found an 
excellent fit with \ 2 — 0.052, which means that our scenario is able to reproduce the lower 
values of the light quark masses as well. The second fit has some qualitative differences in 
comparison with the first one, which reinforces the suspicion that, for given input values 
of the 18 observables, the fit solution in our scenario is not unique. 

Summary: In this paper we have investigated fermion masses and mixings in the 
5*0(10) MSGUT, augmented by a 120-plet of scalars. The main purpose was to show 
that in this setting it is possible to reconcile the type I seesaw mechanism (see Eq. (J6])) 
with a triplet VEV wr equal to the GUT scale of 2 x 10 16 GeV, provided the theory ad- 
mits that the MSSM Higgs doublet Hd is composed mainly of the corresponding doublet 
components in the 126 and 210 scalar irreps — see Eq. (1A11I) ; those are the irreps which 
have no Yukawa couplings. This reconciliation was feasible within the scenario defined 
in points i)-iv), in which we have used symmetries to significantly reduce the number 
of degrees of freedom in the Yukawa couplings — see Eq. ([9]). Within this scenario we 
were able to find an excellent fit for all fermion masses and mixings; in this fit we have a 
hierarchical neutrino mass spectrum^ 

Thus we have obtained the following results for the minimal renormalizable 5*0(10) 
GUT, with Yukawa couplings according to the relation ([T]): 

• It is possible to reproduce the correct neutrino mass scale. 

5 We have also tried fits for the inverted ordering. In that case, the best fit we found has % 2 = 1.8 and 
m 3 ~ 7 x 1CT 6 eV. 
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• Nevertheless, gauge coupling unification is not spoiled. 

• The concrete SO (10) scenario with type I seesaw mechanism, we treated in this 
paper, has 21 parameters, just as the MSGUT with type I+II seesaw mechanism 
and general complex Yukawa couplings. 

It remains to be studied if our scenario allows a sufficient suppression of proton decay. 
In [32] it was shown that the scalar 120-plet plays a crucial role for that purpose; a certain 
texture of the Yukawa coupling matrices — similar to our numerical solution (fTT|) — enables 
that suppression even for large tan/3. 

Acknowledgments: W.G. thanks CS. Aulakh for illuminating discussions and L. Lavoura 
for reading the manuscript. 



A The MSSM Higgs doublets and the mass matrices 

The MSSM contains two Higgs doublets, H d and H u , with hypercharges +1 and — 1, 
respectively. Their corresponding VEVs are denoted by v cos /3 and v sin/3 (v = 174 
GeV), respectively. Neglecting effects of the electroweak scale, these doublets are, by 
assumption, the only scalar zero modes extant at the GUT scale; this requires a minimal 
finetuning condition [TTl [19]. The scalar irreps 10, 126, 126, 210 contain each one 
doublet with the quantum numbers of H d , whereas the 120 contains two such doublets. 
The H d is composed of these doublets [IH] with the corresponding amplitudes [22] ctj 
(j = 1, . . . , 6). The analogous coefficients for H u are denoted by acj. The normalization 
conditions are 

6 6 

E N 2 = EN 2 = L ( A1 ) 

3=1 3=1 

The Dirac mass matrices, taking into account that the 126 and 210 have no Yukawa 
couplings, are given by 

M a = v cos (3 (c a 1 a 1 Y 10 + c a 2 a 2 Y m + (c a 5 a 5 + c a 6 a 6 ) Y 120 ) (a = d,£), (A2) 
M b = v sin f3(c\ ai Y w + c b 2 a 2 Y m + {c b 5 a 5 + c b 6 a 6 ) Y 120 ) (b = u,D), (A3) 

with Yukawa coupling matrices Y%q, Yj^, Y120 and Clebsch-Gordan coefficients deriving 
from the S'O(10)-invariant Yukawa couplings [20], E3] . The absolute values of the Clebsch- 
Gordan coefficients have no physical meaning and some of their phases are convention- 
dependent. With our conventions, the required information reads 

Cd — s~>u — j~JL — (~^-^ 
i — — — t-i , 

n d „u 1 Jt 1 „D 

u 2 — u 2 — 3 u 2 — 3 2 ' 

4 = -eg = 4 = -eg, (A4) 

n d „u 1 J, 1 „D 

°6 — c 6 — 3 c 6 — 3 c 6 ' 

4/4 = Vs. 
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Equations 0A20 and f]A3|) together with this equation lead to the mass matrices fl2])-([5]). 
Furthermore, comparing Eqs. flA2j) and flA3j) with Eq. ffTUl) . we find 



H' = v cos (3 cf\ai\Y w , F' = v cos/3 |a 2 | Y^q, G' = v cos (3 c^a*, + c§olq 
Comparison with Eqs. (|TT i) -f|T4l ) and using Eq. ( 1A4j) delivers the coefficients 



120- 



th = tan (3 
tan /3 



a i 

a 6 - V3a 5 



^6 — \/3q;5 



r> = tan /3 
= 1 - 



«2 

2a 6 



^6 ~~ v^3 a; 5 



, r B = tan /3 



3a 6 + v^as 



(A5) 

(A6) 
(A7) 



Now we want to check the consistency of our numerical solution given by Eq. (TTTT) and 
Table [2j From ^> r u , it follows that 



V3 



a 5 ~ a 6 



cue 



-V3 



«5 



4 tan^ 

Furthermore, using ~ 1, we find the order-of-magnitude relations 

«5 ~ «6 ~ tan/3/ro. 
Then the first of the normalization conditions ( 1A1I) reads approximately 



(A9) 



E 



tan 2 /3 



.2 I- |2 



2 |=r |2 



r H \Oii\ +r F \a 2 \ + 



12 



«6 



>/3 



«5 



+ a 3 + «4 



1. (AlO) 



This means that | ct _y | 2 <C 1 for j = 1, 2, 5, 6. Therefore, the second normalization condition 
is given by 



E 



- |2 

a, 



I - |2 



- |2 



a 3 + a 4 



(All) 



and the brunt of the normalization has to be supplied by the components of H d in the 
126 and 210, which do not couple to the fermions. This is a consistency condition for 
the scenario presented in this paper. 

To translate the condition f)16p into the formalism presented here, we note that \vd\ 2 <C 
|f n | 2 and sin/5 ~ 1 for tan/3 = 10. Therefore, Eq. f)16p effectively checks if the necessary 
condition | a 2 \ < 1 is fulfilled. 

Finally, it remains to see if our numerical solution respects the perturbative regime of 
the Yukawa sector. It suffices to consider the largest elements of the Yukawa couplings 



Y d = —!—-M d , Y u = —^-M u , Y £ = —!—-M e , )) 
v cos [3 vsmp v cos p 



D 



v sin (3 



M, 



D, 



(A12) 



which reside in Y u and Yd. The largest entry in Y u is the 33-element with the main 
contribution from rjjh' 33 /(v sin/3) ~ 0.59. The 23-element with T£,g' 23 /{v sin/3) ~ 0.64 
dominates in Yq. These numbers demonstrate that for our numerical solution the Yukawa 
couplings remain in the perturbative regime. 
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